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I. INTRODUCTION 

The dynamics of almost on-shell heavy quarks Q with mass m much greater than the 
QCD scale Aqcd can be computed in a systematic expansion in terms of several small 
parameters. In the single heavy quark sector, the dynamics is described by heavy quark 
effective theory (HQET), which has an expansion in powers of as{m) and Kqcd/itl [0. 
HQET can be used to compute properties of hadrons such as the B and D mesons containing 
a single 6 or c quark. The dynamics in the quark-antiquark sector is far more complicated 
than in the single quark sector. At low momentum transfer the QQ pair can form non- 
relativistic Coulomb-like bound states, which are the J /ip and T for the cc and hh sectors, 
respectively. It should be possible to describe the dynamics of nonrelativistic heavy quarks 
using a nonrelativistic effective field theory for QCD. A formulation of this effective theory, 
called NRQCD (nonrelativistic QCD), has been proposed by Bodwin, Braaten, and Lepage 
(BBL) 0. The analogous theory for electromagnetism, NRQED, was developed earlier by 
Caswell and Lepage 0. 

Constructing NRQCD has proven to be more difficult than HQET, the complication 
being that there are many scales involved. In HQET, the only two important scales are the 
quark mass m and Aqcd- In NRQCD there two other important scales, mv and mv'^, the 
momentum and energy of the quarks (where v is the typical quark velocity). Momentum 
regions with (energy, momentum) of order {m,m), {mv,mv), {mv'^,mv) and {mv'^ , mv'^) 
are referred to in the literature as hard, soft, potential and ultrasoft, respectively Q]. The 
effective field theory must be able to correctly reproduce phenomena in all of these regions. 

The simplest approach to NRQCD uses a momentum space cutoff to regulate the loop 
integrals. This has the advantage that the physics below the cutoff A is automatically 
correctly taken into account. However, the usefulness of this approach for computations is 
limited since cutoffs break gauge invariance. Furthermore, the theory does not have manifest 
power counting — loop graphs mix powers oi v. If a mass independent subtraction scheme 
such as MS is applied to the BBL Lagrangian, the v expansion breaks down due to unphysical 
poles introduced by the nonrelativistic approximation. There have been many approaches 
advocated to remedy this situation. 

In Ref. 0], it was shown that it was more useful to formulate NRQCD as a theory in 
which ultrasoft modes couple via the multipole expansion. A velocity power counting rule 
for bound states in nonrelativistic effective field theories was formulated in Ref. 0. The 
leading order term in the effective Lagrangian reproduced the form of the propagator in 
the potential regime. To recover the poles in the gluon propagator that correspond to gluon 
radiation, the gluon propagator l/{v'^{k^)'^ — k.^) had to include subleading terms in v, which 
caused problems with the naive velocity power counting rules. In Ref. 0, it was pointed 
out that the usual matching onto NRQCD violated v power counting if the MS scheme was 
used, and it was shown that the problem could be fixed in the single heavy-quark sector 
by using the same matching conditions as for HQET. In Ref. ||], it was demonstrated that 
the multipole expansion is the appropriate generalization of 0] to the two quark sector. In 
Ref. 0, an effective theory was formulated using two different fields for the potential and 
radiation gluons. A problem which arose in this formulation, however, was that it neglected 
soft gluon modes, which are responsible for the running of the coupling below ^ = m. 

In the threshold expansion 0] , the results of NRQCD are obtained directly from QCD by 



expanding graphs about the relevant kinematic regimes (hard, soft, potential and ultrasoft). 
This technique has recently been used to extract the two-loop corrections to top-antitop 
production near threshold with comparative ease p0|-p!3|. However, it is less simple to 



perform renormalization-group improved calculations in this formulation than in a true 
effective field theory (our results in this paper will disagree with the RGE analysis presented 
in |]14|-) The threshold expansion was written as an effective theory by Griesshammer |15[ . 



In the approach advocated by Pineda and Soto ]TB|-P^ the matching onto the effective 
field theory occurs at two stages. Matching between QCD and NRQCD occurs at the scale 
fi = m, while at fi of order the inverse separation between the heavy quarks NRQCD is 
matched onto a new effective theory which the authors call pNRQCD (p for potential). In 
particular, Pineda and Soto argue that the matching between QCD and NRQCD should 
contain only the hard part of loop integrals, and should be performed using HQET Feyn- 
man rules. By performing the matching exactly at threshold, the Coulomb singularity is 
regulated by dimensional regularization, so the one-loop matching conditions are well de- 
fined. Furthermore, the treatment of soft modes is particularly simple in this approach, as 
they just correspond to the running in the theory between m and mv. 

We argue in this paper, however, that the problem with this approach is that HQET 
Feynman rules do not correctly treat the momentum region between m and mv. In partic- 
ular, in |jl8|] it is argued that the anomalous dimension for the electromagnetic current for 
heavy quark production vanishes. While this is true at one loop, at two loops the current 
has a nonzero anomalous dimension [[10| , p!T|| which HQET Feynman rules cannot reproduce. 



In this paper, we construct an effective theory for NRQCD which has a consistent v 
expansion when loops are evaluated in the MS scheme, and which correctly reproduces the 
two loop anomalous dimension of the heavy quark production current. The Lagrangian we 
use is similar to that of [|I^], however, we do not have to introduce as many extra fields (such 
as soft quarks) as in that formulation. Unlike the pNRQCD approach, we argue that the 
correct matching scale onto the effective Lagrangian (similar to that of pNRQCD) is /i = m, 
not n = mv. The added complication which then arises is that soft modes must explicitly 
be taken into account between n = m and fi = mv, in order to obtain the correct running 
of the potential. 

We also introduce a novel renormalization group (VRG) equation in velocity space that is 
used to sum logarithms of v in the effective theory. The VRG represents the invariance of the 
theory under changes in the subtraction velocity u. The formulation of NRQCD presented 
in this paper allows one to include the effects of the running coupling constant in the quark 
potential by using the velocity renormalization group equations, and to simultaneously sum 
soft and ultrasoft logarithms. 

In Sec. 0, we discuss some general aspects of the problem, and in Sec. |III| we introduce 
the fields required in the effective field theory and discuss power counting and loop graphs in 
NRQCD. The VRG is introduced in Sec. |V|, while in Sec. [V| we illustrate the formalism with 
some examples. In particular, we show that we correctly reproduce the two-loop anomalous 
dimension of the heavy quark production current. We defer the complete RGE analysis of 
heavy quark production to a future paper. 



II. p2/2M OR NO p2/2M 

Consider pair production of a QQ pair near threshold by a virtual photon. We are 
interested in the threshold region, where the fermions are nonrelativistic, so that f <^ 1, 
where 



v = jl-— (1) 

is the velocity of the two final state fermions (ignoring for the moment complications from 
confinement effects in QCD). The electromagnetic current in the full theory matches to 

where ^ and x annihilate quarks and antiquarks, respectively. Ignoring for the moment 
non-perturbative effects, there are three relevant scales in the process: the quark mass m, 
the quark three-momentum p = mv, and the quark energy E = mv'^/2. 

An approach to the problem of unraveling these scales was developed in Refs. ||7|,p!6Hl9 



The authors argued that at the scale m, no distinction need be made between energy and 
momentum, since they are both <C m; it is only at the scale mv that they are distinguished 
in the power counting. The correct effective field theory was therefore argued to be identical 
to HQET. The NRQCD and HQET descriptions differ in how they treat these scales. In 
the HQET approach, the kinetic term in C is taken as a perturbation, while in the NRQCD 
approach, the kinetic term is resummed in the propagator. While this violates m power 
counting, it was shown in [^ that as long as the potential is taken to be instantaneous, and 
real radiation is coupled via the multipole expansion, there is a consistent counting in v. 

At one loop, both approaches yield the correct result for the matching onto the external 
current. In the NRQCD approach the Coulomb singularity ~ 1/v in the QQ production 
amplitude is reproduced by nonrelativistic fermions undergoing instantaneous potential ex- 
change, while in the HQET approach the Coulomb singularity is regulated at threshold by 
dimensional regularization. In the latter case the matching condition is given solely by the 
hard part of the graph, obtained by evaluating the full theory at threshold. Both approaches 
give the well-known result for the matching condition, 

ci(/i) = -2Cf. (3) 

At two loops, however, the approaches differ. The two-loop matching onto NRQCD was 
computed by Hoang for QED |jrD|,|r^, and the computation has been recently extended to the 



non-Abelian case by Czarnecki and Melnikov [0], and by Beneke, Signer and Smirnov [|T3 
These authors find 

/I 1 \ m 
C2(/^) = TC^Cp ( ttC'f + T^^A ) In 1- non-logarithmic terms. (4) 

The electromagnetic current has no anomalous dimension in the full theory, which implies 
that in the effective theory, C must have an anomalous dimension at ^ = m, 




FIG. 1. Fermion-antifermion propagation graph. The two create and annihilate a 
fermion-antifermion pair. 



The anomalous dimension is of leading order in the 1/m expansion, and it is straight- 
forward to verify in either A*^ = or Coulomb gauge that the leading order graphs do not 
give a two-loop anomalous dimension in HQET. 

The situation is rather different in NRQCD, which has a v power counting scheme. In 
this case, the 0(1) anomalous dimension arises at two loops due to a 1/f^ enhancement of 
an O(f^) term in the potential: the 1/f ^ Coulomb enhancement is crucial to this result. We 
will compute the anomalous dimension in Sec. [VQ , and show that it correctly reproduces 
Eq. (^. The distinguishing feature between NRQCD and HQET is that HQET does not 
have the Coulomb divergence. By evaluating one-loop graphs exactly at threshold one avoids 
the problem of Coulomb divergences, and this procedure allows one to compute the one-loop 
matching correction Eq. (^ using HQET. However, in two-loop graphs the internal graph 
is not at threshold, and so is sensitive to the Coulomb singularity. The problem in the m 
counting scheme seems to be that unless the p^/2m term is included in the leading order 
propagator, the effective theory cannot correctly reproduce the propagation of a fermion- 
antifermion pair, such as the graph in Fig. |l|, which vanishes in dimensional regularization 
if HQET propagators are used. 

Thus the effective field theory at /i = rra must resum the p'^/2m term in the propagator 
to reproduce the infrared physics of full QCD, and to correctly reproduce the two-loop 
anomalous dimension. Once the p^/2m term is included in the quark propagator, it is also 
necessary to perform a multipole expansion and include a quark-antiquark potential [Q. The 
matching from QCD to an effective theory with potentials is done at the scale fi = m, so the 
potential in the effective theory at fi = m depends on as{m). Since the dominant momenta 
in the static potential are of order mv, one might expect that the relevant coupling is as{fnv), 
and this is borne out by more detailed studies of the quark static potential pO|-p^. One 



therefore requires that the potential generated ai fi = m must run in the effective theory 
below m. This running can be implemented by the inclusion of soft gluon modes, with 
{E,p) ~ {mv,mv), the importance of which was pointed out by Greisshammer [l^. Soft 



gluon modes should be integrated out of NRQCD, since they can never be produced on-shell; 
nevertheless they must be included in the running between m and mv. 

A. Possible Hierarchies 

In addition to the scale mv and m,v'^ the non-perturbative scale Aqcd also plays an 
important role for real quarkonium. Though Aqcd will not play an important role for the 
analysis of this paper, some aspects of its power counting are worth emphasizing. 
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FIG. 2. Plot of Kqqj) / ma1{mv) as a function of m/AqcDi for Uf = 3. 

For very large m, or equivalently, small ^^(tti), one is in the regime Aqcd "C mv'^ -C 
mv -C m, since Aqcd is formally smaller than any power of as- These inequalities are 
only well satisfied for t quarks; for charmonium and bottomium the situation is closer to 
Aqcd ~ f^v"^ or Aqcd ~ ""^"y, and non-perturbative effects become important. Of course, 
the apparent independence of mv, mv"^ and Aqcd for a Coulomb system is illusory. The 
velocity f in a Coulomb bound state is given by solving v = as{mv), 
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where 
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(7) 



which gives v = as{mv) as a function of m/AQCD- In Fig- H, Aqcd/''^«s('^'^) i^ plotted 
as a function of m/AQCD, with v = as{i7iv). Clearly, for large m/AqcB, one can have 
Aqcd <^ ^'^^ <^ fnv <^ m. However, it is not possible to have Aqcd ^ mv^. The maximum 
possible value of Aqcd/^^'W^ is 
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(8) 



(9) 



[Here e = 2.718.] (Of course these values should only be taken as illustrative, since if 
Aqcd ~ fnv^ the system is no longer Coulombic). For the J/\l/, ""^c/Aqcd ~ 3 and for the 
T, rrife /Aqcd ~ 9 so that mf^/AQCD is not very different in the two cases. 

III. THE EFFECTIVE THEORY 



To construct the effective theory, label the total energy E and momentum P of the heavy 
quark by 
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FIG. 3. Momentum space of size mv is divided into boxes of size mv'^. A point in momentum 
space is labeled by p and k. 



P 



E = fc°, 



(10) 



where the three- vector p is of the order of the soft scale mv, and the four- vector k is of the 
order of the ultrasoft scale mv'^. Momentum space of size mv is divided up into boxes of 
size mv'^. The location of each box is labeled by p, and the points within a box are labeled 
by k, as shown in Fig. ^ The variable p is a discrete label, and k is a continuous label. This 
procedure was originally used by Georgi for HQET [^ , where the four-momentum p^ was 
split between mv^ of order m and the residual momentum k^ of order Aqcd? 



pi" = mv'' + /c^. 



(11) 



In HQET, the velocity t; is a discrete label, and fc is a continuous label, so that one sums on 
V and integrates over k p3|. In our case, we will sum over p and integrate over k. 
The quark field iIj{x) in QCD is replaced by 



ip (x) -^ ipp (x) 



(12) 



The label p represents momenta of order the soft scale m,v, and (the Fourier transform of) 
X represents energy and momenta of order the ultrasoft scale mv"^. 

The decomposition Eq. (pH]) is not unique, since one can redefine k ^ k + q, mv -^ 
mv — q, where q is of order k. This redefinition, called reparameterization invariance, 
leads to constraints on the effective field theory, and relates different orders in the 1/m 
expansion [E^. One can make a similar redefinition here. 



k + q, 



p ^ p - q. 



where q is of order mw^. In terms of fields, this transformation is 

^pp (x) -^ e*''-^V'p_q (x) . 



(13) 



(14) 



The application of reparameterization invariance to spinors in HQET was subtle, because of 
the constraint '^ip = ip, that projected out the particle component of the spinor. In our case 
ipp is a two- component spinor whose upper and lower components represent the amplitudes 
to annihilate a quark with spin ±1/2 along a fixed axis. The transformation Eq. (|T^) does 
not affect the spin labels, so the consequences of reparameterization invariance are similar 
to the case of HQET for spin-zero particles. [Spin would enter if the components of tpp 
represented helicity states.] The basic result is that derivatives on ipplx) should be of the 
form zp + V [ID . 

On-shell gauge fields have energy of order their momentum. One can have propagating 
gauge fields with energy and momentum of order mv or of order mv'^, which are referred 
to as soft and ultrasoft gauge fields, respectively. The gauge fields in the full theory are 
replaced by two different fields in the effective theory, momentum-dependent gauge fields, 
Af^{x), and momentum- independent gauge fields A^{x). The fields A!^{x) represent the soft 
degrees of freedom and Af^i^x) represent the ultrasoft degrees of freedom. The total energy 
and momentum of the soft gauge fields is 



P = p + k, E = p" + r, (15) 

and of the ultrasoft gauge fields is k^^, where k is the Fourier transform of the spacetime 
argument x. Note that soft gauge fields are labeled by a four-vector p, whereas quark fields 
are labeled by a three-vector p. Any other light modes (such as light fermions and ghosts) 
in the theory must also be divided into soft and ultrasoft fields, as for the gauge fields. 

The terms in the NRQCD effective Lagrangian describe the interactions of the soft 
gauge fields among themselves, and the interaction of two or more soft gauge fields with 
the fermions. There are no terms that involve the interaction of a fermion with a single 
soft gauge field, i.e. no vertex of the form -ipLA^tjjp, since energy cannot be conserved in the 
interaction. 

The effective Lagrangian for NRQCD can now be written down in terms of the fields ipp 
which annihilate a quark, Xp which annihilate an antiquark, A!^ which annihilate and create 
soft gluons, and A'^ which annihilate and create ultrasoft gluons. The covariant derivative is 
/)M = ^M + igA^" = (L>°, -D), so that D^ = 9° + igA^, D = V - igA, and involves only the 
ultrasoft photon fields. The effective Lagrangian is gauge invariant under ultrasoft gauge 
transformations, those in which the gauge parameter varies on a distance scale l/(mf ^). The 
full gauge invariance of the original Lagrangian is recovered by combining gauge invariance 
in the effective theory with reparameterization invariance. 

The effective Lagrangian in the center of mass frame is 
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FIG. 4. One loop vertex correction in the full theory and effective theory. The momenta are 
related hy P = p + k, P' = p + k' . 

where we have retained the lowest order terms in each sector of the theory. The matrices 
T^ and T^ are the color matrices for the 3 and 3 representations, respectively. The field 
strength tensor igF^'^ = [D^, D"] is constructed only out of ultrasoft gauge fields. If we 
were interested in including the effects of light quarks we would need to add new soft and 
ultrasoft fields for each flavor. The matching at the scale m would then introduce additional 
non-local four fermion terms involving both the heavy and light quark fields. 

In using Eq. (|16|), it is crucial to expand out the term (p — iJD) . The p^ piece is part 
of the leading order Lagrangian that gives the ipp propagator, 



4'"° - ^h 



(17) 



and the terms involving D are treated as a perturbation. This is the momentum space 
equivalent of the multipole expansion written in x space in [||, since the ultrasoft gluons 
do not transfer three momentum to the quarks. This procedure will be justified when the 
velocity power counting rules for the Lagrangian are derived. 

The QCD Lagrangian contains gauge fixing terms and ghost interactions. It is convenient 
to quantize the theory in background field gauge. The background fields can be taken to be 
the ultrasoft modes of the effective theory. The quantum fields represent the quark potential, 
and the soft modes of the effective theory. The effective Lagrangian is gauge invariant with 
respect to the ultrasoft modes, and contains the gauge fixing terms of the original theory 
for the soft modes. One can then gauge fix the ultrasoft gauge fields, to compute loop 
graphs involving the ultrasoft gauge fields. We will use Feynman gauge for both the soft 
and ultrasoft modes. 

The terms in the effective Lagrangian Eq. (|16|) are obtained by matching graphs in QCD 
and the effective theory. The ultrasoft gluon fields A^^ in the effective Lagrangian cannot 
change the momentum label p on the fermion lines, so the single-quark terms in L do not 
change p. The Lagrangian in the single-quark sector is the same as the HQET Lagrangian, 
as pointed out in Ref. [0]. This relation between the NRQCD and HQET Lagrangians holds 
even if loop corrections are included. Consider a loop contribution to a term in the single- 
quark sector of the effective theory, such as that shown in Fig. ^ In the effective theory, 
the incoming and outgoing momenta of the quark, P and P' respectively, are broken up 
into a soft piece p and an ultrasoft piece k, k' . The soft label p must be the same on the 
incoming and outgoing lines, since the ultrasoft gluons do not change p, so the momentum 
transfer in the effective theory is k' — k = P' — P . The matching condition can be computed 
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(a) (b) 

FIG. 5. The Coulomb potential in the full theory, (a), is given by a local two-quark operator 
(b) in the effective theory. 

as the difference between Fig. ^(a) and (b), expanded in a power series in 1/m. The full- 
theory computation is given by computing Fig. ^(a) on-shell, and expanding in powers of the 
external momenta over m. The effective theory contribution is given by evaluating Fig. ^(b) 
on-shell. The on-shell condition in the effective theory is k^ = p^/2m. The intermediate 
fermion propagator 

^ (18) 



^0 + gO _ p2/2m + ie 
is equal to 

1 



g° -|- ie 



(19) 



when the on-shell condition is used.Q The Feynman rules Eq. (|T9|) are precisely those that 
would be used to match from QCD to HQET, and are known not to violate the 1/m power 
counting in the effective theory. Thus the couplings of the ultrasoft gluons in the single- 
quark sector are precisely those in HQET. This was the procedure used to compute the 
HQET and NRQCD Lagrangians at one-loop in Ref 0. 

The Coulomb potential can scatter quark states from one value of p to another. These 
effects are explicitly included in the two-body terms in L. The Coulomb potential is usually 
thought of as a nonlocal two-body operator. However, because of the use of the extra label 
p, the Coulomb potential, 

^^V^tyA (^) ^^ (^) ^t_^ (^) fA^_^ (^) (20) 



is local, and manifestly gauge invariant. The Coulomb potential is obtained as the value of 
Fig. 1^ evaluated in the full theory. The Coulomb potential is proportional to the Casimir 
T^T^, and gives an attractive interaction in the color singlet channel and a repulsive inter- 
action in the color octet channel. 

The leading terms involving the soft gluons are given by matching the Compton scattering 
graphs Figs. ^(a,b,c) in the full theory to the local operator Fig. |^(d) in the effective theory. 
Soft gluons have energy and momenta of order mv, whereas the quarks ipp have energies of 



^Note that we have used the lowest order Lagrangian Eq. ( [T7| ) to derive the propagator, so the 
energy term in the denominator involves the loop momentum q, but the momentum term p^/2r?i 
does not. 
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FIG. 6. The Compton scattering graphs (a,b,c) match on to a local operator (d) in the effective 
theory. Soft gluon modes are denoted by a zigzag line. 

order mv"^ and momenta of order mv. The intermediate quark in Fig. ^(a,b) is off-shell, and 
the Compton scattering graph in the full-theory can be replaced by a local vertex in the 
effective theory, as shown in the figure. The intermediate gluon in Fig. |](c) is also off-shell, 
since energy cannot be transferred from the external gluons to the quark line. Thus the 
interaction in Fig. ^(c) can also be represented by a local vertex in the effective theory. The 
leading order soft interaction vanishes for QED. We will comment on this in Sec. 0. Loop 
graphs involving the soft interaction of Fig. ^(d) are part of the running potential in the 
effective theory. 

The Lagrangian Eq. ( |T^ ) is similar to the pNRQCD Lagrangian constructed by Pineda 
and Soto in Ref . ||2^ , but there are a few important differences. The Lagrangian is given by 



matching to QCD at the scale fx = m, rather than at the scale /i = mv. The Lagrangian also 
contains explicit soft modes. Soft modes are necessary to reproduce the running potential 
in the effective theory. Finally, the Coulomb potential is constructed using a p-dependent, 
but local in x two-body operator, instead of a two-body operator nonlocal in x. The use of 
p-dependent fields is the momentum space analog of the multipole expansion, and simplifies 
the discussion of gauge invariance, particularly in the non-Abelian case. 

A. Power Counting in the Lagrangian 

The NRQCD effective Lagrangian can be used to compute processes to a given order in 
the velocity v. To determine the order in f of a given diagram, it is useful to have a velocity 
power counting scheme, and we will use the one in Table |. This velocity power counting 
scheme differs somewhat from that in BBL, since we have separated the gluon field into 
soft and ultrasoft modes, and the Coulomb interaction. The order in v of the soft gluons 
modes Ap is irrelevant, since they cannot appear as external states in the processes we are 
considering. The final power counting formula for graphs we will derive in Eq. (^O]) holds 
regardless of the order in v chosen for Ap. The lowest dimension operator in the zero-quark 
sector of the Lagrange density is the gauge kinetic energy, which is of order v^. All the 
terms in the field strength tensor F'^'^ are of the same order in v, since D'^ and A^ are both 
of order f ^ . 

The lowest dimension terms in the one-quark sector are 

11 



TABLE I. Velocity counting rules for the effective theory. The electric and magnetic fields are 
those constructed out of the ultrasoft gauge potential Af^. 

p V 

A^ V 

D ^2 

E v^ 

B u^ 



^J(^)|^^°-|^}^p(^) (21) 

which are of order v^. The lowest order Lagrangian Eq. (^TJ) must be used to determine the 
i/jp propagator. Terms that involve the covariant derivative D are of higher order than v^. 
For example, 

^i^^^P (22) 

is of order v^, and must be treated as a perturbation in the effective theory. Each replacement 
of p by D increases the order in v by one. 

The lowest dimension operator in the two-quark sector is the Coulomb interaction. Each 
quark field is of order v^^'^ and l/(p — q)^ is of order l/v"^, so the Coulomb interaction is of 
order a^f^. 

The lowest dimension terms that involve soft gluons are of order OgV^. Consider, for 
example, 

-Ana, Y. i<[4'^;]V^P- (23) 

The two tpp fields are of order v^, the two A^ fields are of order v"^, and 1/q^ is of order 1/v, 
so the vertex is of order a^v^. 



B. Loop Graphs 

The computation of loop graphs using the effective Lagrangian Eq. (|T^) involves some 
subtleties. There are three kinds of loops, which we will refer to as ultrasoft, potential and 
soft, respectively. We will determine the dominant momentum region for each graph by 
studying the pole structure of the diagram, which is what determines the behavior of the 
graph in the MS scheme. 

Consider graphs in the effective theory that involve only a single fermion line, such as 
Fig. §(b). The internal fermion propagators are given by the lowest order term in the one- 
fermion sector, Eq. (|18|). These terms depend only on the soft momentum p, and not on 
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(a) (b) 

FIG. 7. One loop correction to Coulomb Scattering in the full and effective theories. 

the ultrasoft momentum k carried by the gluon hnes. Thus the two fermion propagators in 
Fig. Kb) are 

respectively. These propagators do not depend on the space part of the loop momentum q, 
but they do depend on m. The propagators have poles at energies of order p^/2m or k^. The 
pole positions are set by the external gluon and quark energies, and also by p^/2m, which 
depends on the quark momentum. The fact that the energy poles are determined by the 
momentum (and vice versa), is important, because it relates the soft and ultrasoft scales. 
We will refer to a typical external quark energy as of the order of the ultrasoft scale /^u, and 
a typical external quark momentum as of the order of the soft scale /is- Then loop graphs 
in the one-fermion sector have poles at energies of order /iu and ii^/m. The gauge boson 
propagator is l/[(g°)^ — q^], so the typical momenta in the loop are of the same order as the 
gluon energy, q ~ q°. Loop graphs such as Fig. ^b) in which the ultrasoft momentum q 
is integrated over will be referred to as ultrasoft loops, and when evaluated in dimensional 
regularization, are dominated by energy and momenta of order /xu or fis/rn. 

The use of p^/2m rather than (p + q) /(2m) in the propagators removes the problem 
of the breakdown of the effective theory due to poles of order m in loop graphs. The power 
counting scheme of Table | in which p is of order mv, but q is of order mv'^ requires that 
(P + q) be expanded as p^ + 2p ■ q + q^, with p^ included in the fermion propagator, and 
the higher order terms 2p • q + q^ treated as vertex insertions. 

Consider a graph that involves a potential loop, such as the one-loop correction to 
Coulomb scattering. The graph is shown as Fig. |^(a) in the full theory, and as Fig. |^(b) 
in the effective theory. The external fermions have energy E and momenta ±p, p'. In the 
effective theory, the external fermions are labeled by the soft momentum ±p, p', and the 
ultrasoft momentum {E,0). The intermediate fermions in the effective theory have soft 
momentum ±g, and ultrasoft momentum {E, 0) ± k. The graph in the effective theory is 

°^^y^^(p-q)2 (p'-q)2 A;0 + E-qV2m -A;0 + E-qV2m' ^^^^ 

There is an integral over the ultrasoft energy fc°, the ultrasoft momentum k, and a sum 
over the soft momentum q. Recall the decomposition of momentum space shown in Fig. |^. 
Summing over q and integrating over k is equivalent to integrating over the entire momentum 
space. Thus one can replace Eq. ( P3[ ) by 
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FIG. 8. One loop graph involving the soft gauge fields. This graph contributes to the running 
potential in the effective theory. 

4 1 1 1 1 

^(P - q)^ (P' -q)^q^ + E- qV2m -q^ + E - qV2m' ^^^^ 



where 



dk° ^dq°, T. fd^-^ f d<i (27) 



The loop graph is dominated by g° ~ E, and q ~ \rrnE, which are ultrasoft, and soft, 
respectively. This is consistent with the picture that g" ranges over a box of size mv^ and q 
over a box of size mv. It again shows that one cannot treat the soft and ultrasoft scales as 
independent of each other; they are related by (soft) ~ ultrasoft x m. 

Finally, consider a graph such as Fig. ^ that involves a soft loop. It gives a contribution 
of the form 



OC y fd% —-:r —— :f 2 :j , (28) 

Vy D-d' 2 D-D' 2 (nO?-a^ r„o^2_/ , a2' V / 



where we have used the last interaction in Eq. (|T6|) for the soft vertices. The sum on q is 
over a four-vector. As for potential loops, one can make the replacement 

yfd'k^fd'q, (29) 

where the replacement must be done for all four components of q, since soft gluons carry a 
four-vector label q. It is straightforward to see that Eq. (pS]) is dominated by q^ and q of 
order mv, which is consistent with using the replacement Eq. (^) for all four components 
of q. 



C. Power Counting Formula for Loop Graphs 

It is now straightforward to derive a power counting rule for an arbitrary graph in 
NRQCD. A given graph has L^ ultrasoft loops, Lp potential loops, and Ls soft loops. 
These can be determined from the structure of the diagram in a systematic way. The total 
number of loops is Lu + Lp + Ls- Now delete all the ultrasoft lines from the graph. The 
remaining number of loops is Lp + Ls ■ Finally, delete all quark lines from the diagram. The 
remaining number of loops is Ls- An example of loop counting is shown in Fig. ^ 

Let Vk denote the number of vertices of order v'^ in a given graph. For example, ip^iD'^ijjp 
is a vertex of type V5 , since it is a one-fermion vertex of order v^ . It is convenient to break 
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^^,jSf(rJW«-«r^^^ 




(a) 




(b) 



(c) 

FIG. 9. An example of counting ultrasoft, potential and soft loops and vertices. Graph (a) 
has 9 vertices and 6 loops. Deleting all ultrasoft lines gives graph (b), which has 9 vertices and 
4 loops. Deleting all fermion lines from graph (b) gives graph (c), with 4 vertices and 2 loops, 
which is the number of soft vertices and loops. There are 5 potential vertices and 2 potential loops 
(the difference between (b) and (c)), and 4 ultrasoft vertices and 2 ultrasoft loops (the difference 
between (a) and (b)). The graph (c) has two connected soft components, Ns = 2. 
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up the vertices into ultrasoft, V^ , potential V^ and soft V^ . The ultrasoft vertices V^ 

(p) 
are those vertices in Vk that involve only ultrasoft fields, the potential vertices Vj^ are those 

(S) 

with at least one fermion field and no soft fields. The soft vertices V^ involve at least one 
soft field. An example of vertex counting is shown in Fig. ^ 
A diagram is of order v^, where 



k 



k [vi^^ + VP + Vi^'^] - 5If - AIs - 8Iu + 8Lu + 5Lp + ALs- (30) 



The first term simply adds up the powers in v of all the vertices. Each internal quark line 
eliminates two ip-p fields, and gives a factor of the fermion propagator 1/{E — p^/2m), which 
give a net factor of l/(f^f ^) = v~^ . This gives the term —5If, where Ip is the number of 
internal fermion lines. Each internal soft line eliminates two Ap fields, and gives a factor of 
the gauge propagator l/((p°)^ — p^), which gives a net factor of l/(f ^f ^) = f ""*. This gives 
the term — 4/^, where Is is the number of internal soft lines. Each internal ultrasoft line 
eliminates two A^ fields, and gives a factor of the gauge propagator l/{{k^y — k^), which 
gives a net factor of l/(f^f^) = v~^. This gives the term —8Iu, where /[/ is the number of 
internal ultrasoft lines. 

Ultrasoft loops are dominated by energy and momentum of order mv'^, and give give a 
factor of f ^ for each integral, so that one gets a factor of v^ for each loop. Potential loops 
are dominated by energy of order mv'^ and momentum of order mv, and give a factor of 
f ^ for each time integration, and v for each space integration, for a net factor of v^ per 
loop. Soft loops are dominated by energy and momentum of order mv, and give a factor 
of V for each integration, for a net factor of v^ for each loop. These contributions give the 
8L(7 + 5Lp + 4Ls term in Eq. 

The identity 



E [vf ' + vr + vr 

k,n 



Is-Iu + Lu + Lp + Ls = 1, (31) 



is the usual relation that the Euler character of a connected graph is unity. An analogous 
relation holds for the graph with all ultrasoft lines removed. We will assume that the graph 
remains connected when ultrasoft lines are removed, which is true for any process in which 
momentum of order mv is transferred between the two fermion lines. The relation for the 
graph with ultrasoft lines removed is 



E [v^^ 



(P) , T/('5) 



+ V^''> -If-Is + Lp + Ls = 1. (32) 



For Ip to be equal in the two relations Eq. pT| ) and (0), it is important that one not erase 
vertices where the gluons couple to the fermions (see Fig. P), so the total number of vertices 

(p\ fQ\ 

is given by V^ + V^ . Finally, one has the Euler character relation for the graph with all 
ultrasoft and fermion lines removed, 



Y.V^'^-Is + Ls = Ns, (33) 

k 

where Ns is the number of connected components in the soft graph. 
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Eliminating lu. If and Ls between Eqs. (|30|-p3D gives the result 



(5 = 5 + E [(A: - mP + ik- 5)VP + {k - 4)V, 



(S) 



-N.. 



(34) 



A given soft vertex in the Lagrange density has the generic form 

1\^ 



ii^priA,)' (Ar p' (DY [-)' 



(35) 



where ijjp represents any quark or antiquark fields or their conjugates. The term Eq. (^ 
has dimension four, 



and is of order f '^, where 



4= -a + b + c + d + e- f, 



(36) 



k = -a + b + 2c + d + 2e. 



Subtracting these two relations gives 



fc-4 = c + e + /. 



(37) 



(38) 



The vertex can only have positive powers of D, A^^, and 1/m, so c + e + / > 0. [Note that d 
need not be positive.] For a soft vertex, let a = c + e + / > 0. Label the soft vertices V^ 

(S) 

by the values of b and a, so that they are denoted by V^^ , where k = A + a. Then one finds 



r{S) 



E(^-4K^ = E<. 



6,0- 



Substituting this result into Eq. (|34D gives the power counting formula 
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5 + E [(^ - 8)V^r + (k - 5)VP] + E ^C 



(s) 



Ns 



(39) 



(40) 



This is an important result. All terms in the zero-fermion sector have k > 8, and in the 



(u) 



r{S) 



nonzero fermion sector have a > 0. Thus the contributions of {k — 8)V^ , and cV^^ are 
each positive. There can be negative powers of v from both soft and potential exchange, 
but these come with compensating factors of as- The Coulomb interaction is fc = 4 vertex 
in V^^\ but is of order a^. Iterating the Coulomb interaction n times produces terms of 
order (a^/w)", so for v ~ a^, the Coulomb interaction must be summed to all orders. The 
Ns term is negative, but each soft component must contain at least one power of a.^- 
As an example of Eq. 



(|3^), consider the electron self-energy diagram Fig. |T0| which 
produces the Lamb shift. The graph contains two p ■ A vertices, each of which is of order v^, 
and one ultrasoft loop, so the net power is 5 = 5 + (6 — 5) x 2 = 7. The graph has a factor 
of as from the two gauge couplings, and so is of order a^f ^ compared to the leading term 
in the single-fermion sector (which is of order v^), i.e. it is of relative order a^f^. Similarly 
the soft loop graph in Fig. ^ has a = 0, and Ns = 1, and so is of order a'^v'^. This is the 
same order in v, but one higher order in as than the Coulomb interaction. 
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FIG. 10. Graph contributing to the Lamb shift. The dot represents a p • A interaction. 
IV. VELOCITY RENORMALIZATION GROUP EQUATION 



Loop diagrams in the effective theory can be divergent, and the effective theory is renor- 
mahzed using the MS scheme. A scale parameter fi\j is introduced when one analytically 
continues the Lagrangian from four to D = 4 — 2e dimensions, as in conventional dimensional 
regularization of gauge theories. There is a subtlety in the use of dimensional regulariza- 
tion in NRQCD. In evaluating loop graphs with potential loops, we made the replacement 
Eq. (p?!). In D dimensions, the relation should read 



de-^dq\ E/^^'^'k 




i-D 



d 



D-l 



q- 



(41) 



\i~D 



The factor of [jio,/ iJi\]Y~^ is needed to ensure the correct dimensionality on the two sides 
of Eq. (0). The integral over (ik is over a volume of the order of /Xy"^, since the typical 
range of integration is of order /xu- The integral over t/q is over a volume of the order of 
/Xg ~^, since the range of integration is of order /ig. The number of terms in the sum on q on 
the left-hand side of Eq. ( |1T| ) is the ratio of the two volumes in four dimensions, (/is//^u) ■ 
Away from four dimensions, this number does not properly account for the momentum space 
volumes on the two sides of Eq. (^l|) , and the additional factor of (/is//Uu) ~ is needed. The 
mismatch of dimensions in D 7^ 4 occurs only for the space part of the integral. The factor 
above is correct when the NRQCD integrals are done in the conventional way, by first doing 
the k^ integral using the method of residues disregarding the contour at infinity, followed by 
the k integral in 3 — 2e dimensions. Similarly, for soft loops Eq. (^) should be replaced by 



ajd^'k 




4-D 



/.-. 



(42) 



The effective theory renormalized in the MS scheme has two /i parameters, /is and /iu- 
However, the two parameters are not independent, since the soft and ultrasoft scales are 
related, /ig = r/i/iu- It is better to think of the parameters as /ig = mv and /iu = mz/^, where 
v is the subtraction point velocity. One can now derive a new kind of renormalization group 
equation for the effective theory, since the bare theory is independent of the subtraction 
velocity v. This velocity renormalization group equation can be used to scale the coefficients 
in the effective theory from the matching scale /ig = m, /iu = "^ to /ig = mv, /iu = mw ^, i.e. 
from V = 1 to V = V. 

The velocity renormalization group equation addresses an important point about the 
effective theory, the simultaneous existence of two related scales /ig and /iu- Loop graphs 



involving gluon loops will typically have logarithms of the form In fi/E, where E ~ mv'^ 
is the typical photon energy. The scale fi is equal to /iu, since fj,\j is the scale parameter 
introduced in D dimensions. Loop graphs containing four-fermion terms, (the potential 
and soft loop graphs discussed earlier), typically have logarithms of the form In fi/ymE or 
In/i/p. These graphs use the replacement Eq. (pj), and so have have a factor of 




= /^r (43) 

for each loop, where the first factor is from Eq. (^TJ), and the second factor is the conven- 
tional factor for each loop in dimensional regularization. Thus in potential and soft loops, 
logarithms are of the form In^s/V^nE or In/is/p- The radiation and potential logarithms 
are 

ln^ = ln^, ln^ = ln^. (44) 

The choice of renormalization point u = v ensures that both logarithms are simultane- 
ously small. Thus using the velocity renormalization group equation from i/ = 1 to i/ = v 
simultaneously sums the logarithms involving the soft and ultrasoft scales. The velocity 
renormalization group equation is the MS equivalent of using an energy cutoff mu'^ and a 
momentum cutoff mi> in a hard cutoff scheme. The VRG allows one to have a static poten- 
tial with an effective coupling constant asimu), and radiation corrections with an effective 
coupling constant as{'mi>'^). 

In a conventional renormalization group approach, one would scale the effective theory 
from fi = m to fi = mv, and then down to /x = mv'^. The VRG differs from this in 
an important way, because it uses a subtraction velocity rather than a subtraction scale. 
Scaling the theory from z/ = 1 to z/ = f is equivalent to simultaneously scaling potential and 
soft graphs from m to mv, and radiation graphs from m to mv^. The scale mv and mf^ 
are coupled in the theory, and this coupling of scales is better treated using a subtraction 
velocity rather than a subtraction scale. 



V. EXAMPLES 

The formalism we have developed will be applied to three illustrative examples in this 
section, the one-loop correction to the static potential, integrating out a heavy fermion, and 



the two-loop anomalous dimension of the production current [10, 11 



A. Box Graph and the Static Potential 

The first example we consider is the renormalization of the static potential at one-loop. 
At tree-level, the fermion-fermion scattering amplitude is reproduced in the effective theory 
by a local operator, the Coulomb vertex in Eq. (^). At one- loop, the QCD diagrams that 
contribute are shown in Table |I| above the horizontal line, and the graphs in the effective 
theory are shown below the horizontal line. 
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TABLE II. One-loop correction to quark-antiquark scattering. The color factors listed in the 
table are for the color singlet channel. The column labeled QCD gives the contribution of the 
graph evaluated in QCD in the MS scheme in Feynman gauge, from Titard and Yndurain p6|,|27|. 
Only the logarithmic contributions are given. A is a gluon mass used as an infrared regulator, k 
is the momentum transfer, and fi is the scale parameter of dimensional regularization. The gluon 
vacuum polarization graph includes the contribution from the ghost loop. The column labeled 
HQET gives the contributions of the diagrams in HQET, where m — > oo. The breakup of the total 
diagram into soft and ultrasoft contributions is given in the last two columns. The two diagrams 
below the horizontal line are the contributions to quark-antiquark scattering in the effective theory. 
The complete set of ultrasoft diagrams is show in Fig. Illl 



Diagram Color Factor QCD HQET Soft Ultrasoft 



„/%. 



J2l 






c. 



(Cp — gC'A 



Cf 



[Cf — gC'A 



C, 



Ca 



21nt 



21nt 



-21n^ -21n^ 

k k 



In ' 



-2 In. 



In ^ 2 In ^ 



fln^ 

2 m 



-|ln- -|ln- 



-2In^ 



2 In 



2In^ 



-2In- 



-2In- 



2In^ 



-|ln^ 




Ca 



liln^ 

/i 
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The sum of the QCD diagrams gives the net one-loop contribution to the static potential 
(in the color singlet channel), 



V(k) = _l!L^i(/^ 

|k| 



6 TT H 



(45) 



where only the logarithmic term has been retained. This gives the well-known result that 
the potential can be rewritten as 

r(k) = '""-f'^ r.T, (46) 

|k| 

where T ■ T = —Cp in the color singlet channel. 

The box graph (the first diagram in Table |Tp was computed using the threshold expan- 
sion by Beneke and Smirnov ||^, and it is interesting to see how the effective field theory 
reproduces the various contributions. The hard part of the box graph is the matching condi- 
tion between QCD and NRQCD, i.e. a local two-fermion operator in the effective Lagrangian, 
and is also equal to the difference between the graphs computed in QCD and HQET. The 
potential part of the box graph is reproduced in the effective theory by the one-loop contri- 
bution from the iteration of two Coulomb interactions Fig. |^. The HQET value of the box 
diagram, 2Cf In A//c, is both infrared and ultraviolet finite if a gluon mass A is introduced as 
an infrared regulator. This total contribution is split, in the the threshold expansion, into 
an infrared divergent soft contribution 2CFlnfi/k, and an ultraviolet divergent ultrasoft 
contribution 2Cf In X/fi. The sum of the soft and ultrasoft contributions has no 1/e pole, 
since the ultraviolet and infrared divergences cancel .0 

The ultrasoft contribution vanishes in the full and effective theories, if the infrared di- 
vergence is regulated by dimensional regularization. The non-trivial contribution is the soft 
part of the box graph. The importance of this contribution was pointed out by Griessham- 
|T5| , who argued that one needed to introduce soft gauge fields, as well as soft quark 



mer 



fields. In our approach, the soft part of the box graph is reproduced in the effective theory 
by the loop graph shown in Fig. p, where the interaction vertex is from the last two lines of 
Eq. (|TBp. It is not necessary to introduce soft quark fields, as advocated by Griesshammer. 

In QED, the soft part of the box and crossed-box graph cancel. This is consistent with 
the effective Lagrangian Eq. (|16D, where the soft interaction vertex vanishes for QED. In 
QED, the soft modes can be integrated out directly at u = 1, and replaced by local operators 
at the scale ^ = m. This approach does not resum the logarithms off in the potential (which 
are absent for QED), and so is not a satisfactory procedure for QCD. 

The effective theory correctly reproduces the soft and ultrasoft contributions to the static 
potential. The soft vertex in Fig. ^ is computed from the Compton scattering graphs in 
Fig. p. Fig. 1^ reproduces the sum of the soft part of the box, vertex and vacuum polarization 



^The cancellation is not really between an ultraviolet and infrared divergence. In the effective 
theory, there are tadpole graphs which are zero in dimensional regularization, and have the form 
of a difference 1/e — 1/e between an ultraviolet and infrared divergence. One cannot characterize 
a 1/e pole as an infrared or ultraviolet divergence if tadpole graphs are set to zero. 
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FIG. 11. Renormalization of the four-quark operator representing the quark- antiquark potential 
by ultrasoft gluons. 



corrections to fermion-fermion scattering. We have seen that the soft vertex is proportional 
to the commutator of two gauge fields, so the soft graph Fig. | is proportional to Ca- This 
automatically implements the cancellation between the various Cp contributions to quark- 
antiquark scattering in the full QCD calculation. An explicit computation of Fig. ^ gives 
the contribution 



47ra.T ■ T 



a. 



■K 



Ca 



. 5\ k 
6/ /^s 



(47) 



to the scattering potential, where the first term (1) is an infrared divergent contribution, 
and the second (5/6) is an ultraviolet divergent contribution. Note that the soft graph 
is infrared divergent even if a gluon mass is used as an infrared regulator. The infrared 
divergent contribution of Eq. (^Tf) is converted to an ultraviolet divergent contribution if 



tadpole graphs are included. Equation ( |77| ) agrees with the sum of the soft contributions 
listed in Table ||. 

The ultrasoft contributions in the QCD theory add up to zero. The ultrasoft contribu- 
tions in the effective theory are the renormalization of the local four-fermion quark- antiquark 
potential, shown in Fig. ^. Each graph is ultraviolet divergent, and proportional to lnfi\j/k. 
The sum of all the graphs is zero, in both the singlet and octet channels. In the effective 
theory, ultrasoft radiative corrections do not renormalize the quark-antiquark potential. 
They do, however, cause mixing between the leading order potential, and corrections to the 
potential suppressed by powers of 1/m. 

The quark-antiquark potential is VRG invariant. This implies that the coupling a^ of 
the soft gluons must satisfy a /3-function equation, where the /3-function for the VRG is the 
same as the conventional one, since /is = mz/, so that fisd/dfis = vd/dv. In other words, the 
quark-antiquark potential takes the form Eq. (|^), with the coupling constant renormalized 
at the soft scale /i = /is = rnv. Choosing mv = |k| sums the leading logarithms, and gives 
Eq. (IID. 

The effective theory has performed an interesting rearrangement of the terms in the ra- 
diative correction to the static potential, compared with those in the corresponding HQET 
computation in Feynman gauge. In the HQET computation, the vertex and wavefunction 
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FIG. 12. A light fermion bubble contributions to an ultrasoft loop. 

renormalization graphs contribute —Ca In A;//i, and the vacuum polarization graph con- 
tributes an additional — (5/6)Cyilnfc//i. The box diagram is finite, and does not contribute 
to the running of the static potential. The vertex and wavefunction graphs involve ultrasoft 
loops, whereas the vacuum polarization graph involves soft loops. In the effective theory, 
the box graph contribution In k/\ is broken up into a soft piece. In k/ ^, and an ultrasoft 
piece, In/i/A. The ultrasoft piece cancels the vertex and wavefunction graphs, leaving the 
soft contribution, so that the static potential depends only on /is, and not on /iu-Q The 
breakup of the box graph into soft and ultrasoft can be done without double counting in a 
mass- independent scheme such as MS @]. 

Terms in the one-fermion sector of the theory are renormalized by ultrasoft gluons. The 
ultrasoft gluon coupling constant is renormalized due to their self-interactions. All these 
graphs can be computed as for HQET ||^,^,^. The VRG anomalous dimension is twice 
the usual anomalous dimension, since vd/dv = 2fi-[jd/dfi\j. 

There are relations between the soft and ultrasoft gluon couplings in the effective theory. 
The two couplings have a /5-function that is related to the QCD /5-function, so that the soft 
coupling is Osoftl'^) = asimu) and the ultrasoft coupling is Ouitrasoftl'^) = «s("^'^^), where 
as{^) obeys the usual renormalization group equation. This can be verified by explicit 
computation in the effective theory. 

B. Integrating out a light fermion 

Consider the NRQCD theory with an additional fermion \E' of mass m<i,, with Aqcd *^ 
rriqi <^ m. At the scale m, one can match from QCD to an effective theory that contains, in 
addition to the fields we have been discussing, the fermion \E'. At the scale m, the fermion 
\E' behaves like a massless particle, so the effective theory contains soft and ultrasoft fermion 
fields for \E', "^qix) and "^{x), respectively. The VRG equation is used to scale the theory 



below m. At the velocity mu'^ = m^,, i.e. u = Jm^/m, the ultrasoft fermion modes \E'(a;) 
can be integrated out of ultrasoft loops such as Fig. IT2|, and at the velocity mi> = mq,, i.e. 



v = nixst/m, the soft fermion modes '^q{x) can be integrated out of soft loops such as Fig. |T3 
This sums logarithms of m/m^ in both soft and ultrasoft loops. 



^This is true even when the running soft gluon coupling as{mi^) and ultrasoft gluon coupling 
as{mi>'^) are used. 
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FIG. 13. A light fermion bubble contributions to a soft loop. 

C. T^vo-loop running of the production current 

A highly non-trivial check of the effective theory is the computation of the two-loop 
running of the production current. This was computed by Hoang for QED |iy, |1^, and 



the computation has been recently extended to the non-Abelian case by Czarnecki and 
Melnikov [|I^, and by Beneke, Signer and Smirnov ||13|| . 

Consider production of a QQ pair near threshold by a virtual photon. The electro- 
magnetic current in the full theory matches to the effective current Eq. (H) in the effective 



theory, where the two-loop matching condition |T^-|T3| between the full and effective theo- 
ries is given in Eq. (^). The electromagnetic current has no anomalous dimension in the full 
theory, which implies that in the effective theory, C2 must have an anomalous dimension. 



at /i = m. 

The anomalous dimension Eq. ( ^8] ) is of order alv^ , so we need to compute all diagrams 
of this order in the effective theory. The interactions needed are in the quark-antiquark 
potential to order f ^ in the color singlet channel. The potential in the center of mass frame 
for the process Q(p) + 0(~P) "^ Q{v') + 0(~P') with momentum transfer k = p — p' is 
(borrowing the notation of Titard and Yndurain [^,|2j 



TTilKl rn^ rn^ m? m? 



where 



1 


+ 


1 2 

|P 


k^ 


m2|k|^ 






S 






A(k) 



^' . V2 , Vhf^^ , Vls ,,,. , Vt^ 



2 ' 
.S-(kxp) 



(50) 

(51) 



k2 ' 

T = o-i-(T2-^(k-o-i)(k-o-2) (52) 



and the coefficients we need are 



Vc'(/i = m) = —4:7ras{m)CF, 

1 

V2{n = m) = 0, 



V|k|(/i = m) = al{m)CF ( -Cp - Ca 



'■> 



Vft/(/i = m) = , (53) 
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(a) (b) 





(c) (d) 




(e) 

FIG. 14. Diagrams which contribute to the two- loop anomalous dimension of the production 
current. The dot and box represent terms in the tree-level potential at order 1 (Vc), and v'^ {Vc, 
V2, Vhf, Vls, Vt) respectively. The (81 is an insertion of the order t;^ term proportional to a'^/m |k| 
in the one- loop potential (Vjki). The x is an insertion of the p^/SmP' kinetic energy correction. 

using tree-level matching at /i = m for Vc^hf and one- loop matching at /i = m for V\\f^\. 
Note that the 1/ |k| and |p| / |k| terms have the same coefficient, which follows from repa- 
rameterization invariance. In addition to the above terms, we also need the p'^/^m? kinetic 
energy correction in the Lagrangian, whose coefficient is also fixed by reparameterization 
invariance. 

The diagrams which contribute to the two-loop anomalous dimension of the production 



current in the effective theory are shown in Fig. IT^ The contribution from the diagrams to 



the anomalous dimension will be called 7a-7e, respectively, where 



dC dC 
/is-i— = V—- = 7. (54) 

d/is du 



One finds that 



7a = 0, 

le = y V^ki (55) 

and the total anomalous dimension is 
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^ ^ ^ ^ ^^S Vc{V2 + 2Vhf) , TT^l^i,! /.«N 

7 = 7a + 7fe + 7c + 7d + 7e = — ^ ^^ - + ^^ (56) 

which agrees with the known result for QCD a.t fi = m, Eq. (^), when the coefficients 
Eq. ( p3D are used. 

To complete the renormalization group analysis of the two-loop anomalous dimension, 
one needs the running values for Vc, Vhf and V\k\, which will be presented elsewhere. Note, 



however, that our analysis disagrees with the results of [Q, in which the running of these 
terms in the potential was neglected. 

VI. CONCLUSIONS 

We have discussed a formulation of nonrelativistic QCD that can be used with a mass- 
independent subtraction scheme such as MS. The effective theory passes several non-trivial 
checks: (1) it has a consistent v power counting expansion, (2) it correctly reproduces the 
running of the quark-antiquark potential at one loop, and (3) it correctly reproduces the two- 
loop running of the production current. The effective theory allows one to simultaneously 
treat the momentum regions of order mv and mv'^. 

In the way we have formulated the effective theory, all large logarithms have been summed 
at the scale u = v. At this point, one can integrate out the soft and modes, and at the 
same time switch from a theory of quarks and antiquarks to a theory of quarkonia, i.e. to 
an effective Lagrangian representing quarkonia interacting with background color fields, as 
ffist studied by Voloshin |30] and Leutwyler [^ . 



The methods described here should also be applicable to other nonrelativistic field the- 



ories, such as those describing nucleon-nucleon scattering at low energies [p^ 
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